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We compute the amplitude of deeply virtual Compton scattering (DVCS) using the calculus 
of QCD string operators in coordinate representation. To restore the electromagnetic gauge 
invariance (transversality) of the twist-2 amplitude we include the operators of twist-3 which 
appear as total derivatives of twist-2 operators. Our results are equivalent to a Wandzura-Wilczek 
approximation for twist-3 skewed parton distributions. We find that, for the tensor amplitude 
Tf_tv , this approximation gives a finite result for the longitudinal polarization of the virtual photon 
and a divergent result in case of transverse polarization. However, the divergent part has zero 
projection onto the polarization vector of the final real photon. 
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I. INTRODUCTION 



The deeply virtual Compton scattering (DVCS) process, in which a highly virtual photon 7* (91) produces a real 
photon 7(^2 = Qi + at small invariant momentum transfer t — r'^, is receiving a lot of attention as a potential 
source of new information about the nucleon structure in terms of so-called "skewed" parton distributions (SPD's) 
The first experimental observation of DVCS has been already reported by the ZEUS collaboration ||]. 
Originally, only the purely twist-2 contribution to the tensor amplitude T^^ amplitude was calculated , which 
does not include subleading terms proportional to the transverse part rj^ = A of the momentum transfer. As a 
result, the twist-2 amplitude formally violates the electromagnetic (EM) gauge invariance in terms linear in A. 
Guichon and Vanderhaeghen proposed to add an ad hoc 0(A) term producing an "improved" expression which 
is EM gauge invariant up to O(A^) terms. Their prescription was recently supported by several groups which 
derived this term in a regular way as a kinematical twist-3 contribution. Furthermore, Anikin et al. Q, using the 
momentum-space collinear expansion, have obtained expressions which include (in case of the pion target) all the 
relevant twist ~3 operators. The DVCS amplitude to accuracy Xj \J —q^ was also calculated by Penttinen et al. ^ in 
a parton model approach. Within the light-cone expansion framework, Belitsky and Miiller |^ analyzed both quark 
and quark-gluon contributions and demonstrated that, to get a gauge invariant result up to terms of order t/qf, it is 
sufficient to retain only the part of the twist-3 SPD's which is obtained by Wandzura-Wilczek (WW) type formulas 
from the twist-2 distributions. Recently, Kivel et al. established that the WW-approximation expression for 
the transverse polarization of the virtual photon diverges. The mathematical aspects of twist decomposition were 
discussed by Bliimlein and Robaschik , Geyer and Lazar . 

Our goal is to analyze the DVCS amplitude within the QCD string operator approach of Balitsky and Braun p3| , 
which proved to be a powerful tool to investigate the higher-twist effects. Here, we consider only the kinematical 
twist-3 terms. In this sense, our results are equivalent to the WW approximation. In addition to offering an 
alternative derivation of this approximation, we incorporate the formalism of double distributions [^,0 which 
provides a simple way of deriving relations between new SPD's. For instance, the fact that the WW approximation 
for the amplitude T^^ gives finite results for longitudinally polarized photon, but diverges in the case of transverse 
polarization, can be easily understood on the basis of the formulas relating these SPD's and the basic twist-2 DD's. 
We also note that the divergent part of T^i^ has a vanishing projection onto the polarization of the final real photon. 



II. DVCS AMPLITUDE 



Generalities. The virtual Compton scattering amplitude is derived from the correlation function 

T^.=^J d^x J d% e-^(«^^)+*(«-^\p2|T{J^(x)J.(2/)} \p,) , (1) 

where J'^{x) is the electromagnetic current operator. Due to current conservation c)^Jp(x) = 0, this function is 
transverse with respect to the incoming and outgoing photon momenta: T^^ = , (72 ^a"' = . It is convenient to 
switch to symmetric variables q = (qi + q2) /2, r ^ q2 — qi and p = {pi +p2)/2. Then, the transversality conditions 
convert into two relations 

connecting the symmetric Ts^^^^ = {T^^i, + Tu^)/2 and antisymmetric T^^^^ = {T^i, — Ti,^)/2 parts of T^'^ . In the 
r = forward limit, the two relations decouple to give the DIS transversality conditions q^Ts^^^^^^ = 0, q^T[^^] = 0. 

In DVCS, the initial photon is in the Bjorken kinematics {— — >■ 00, (piqi) 00, xb = — /[2(pigi)] fixed} and 
the final one is real g| = 0. Since <z| = + 2(gir) + 1, the momentum transfer r in this process should have a large 
component in the direction of p, with (rqi) close to — (7i/2 for small t. The size of this component is characterized 
by the skewedness parameter rj = {rq)/2{pq). For DVCS, rj coincides, up to 0{t) terms, with the generalized 
Bjorken variable ^ = —q^/2{pq) . Hence, the momentum transfer may be split into the component parallel to p and 
the remainder A 
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r = 2^p + A, 



(3) 



which in the t = 0, = limit is transverse both to p and q: (Aq) = —t/i, (Ap) = —2^p^. 

Coordinate representation. An efficient way to study the behavior of Compton amplitudes in the Bjorken limit is 
to use the light-cone expansion for the product Ilf^^{x,y) = iTJf^{x)J^{y) of two vector currents in the coordinate 
representation. Following Balitsky and Braun we start from the formal light-cone expansion in terms of QCD 
string operators (with gauge links along the straight line between the fields which, for brevity, we do not write 
explicitly). The leading light-cone singularity is contained in the "handbag" contribution 

n^,(z I X) = ^l^s^p,,0,{z I X) - e^p,„05a{z I X)| , (4) 
where Sf_,p^a = gp-pQi^a - g^L^Qpa + Qtiagyp, X ^ {x + y)/2, z = y - x, and 

= ^ [^{X^z/2)^,i:{X + z/2) - {z^-z)] , 
05a{z \X)^^ [i'{X - z/2)7,75V'(^ + z/2) + {z ^ -z)] . (5) 

Twist-2 part. The string operators in Eq. (Q) do not have a definite twist. The twist-2 part is defined by 
formally Taylor-expanding the string operators in Eq. in the relative coordinate z and retaining only the 
totally symmetric traceless parts of the coefficients in the expansion: 

[^{X - z/2)7,V(^ + z/2)]*""*"' = E ■ • ■ ^i^) [fW Dp,... B^„} -traces] i,(X), (6) 

and similarly for the operator with 7^75. As shown in Ref. [p^ , "symmetrization" and "subtraction of traces" can 
be carried out directly at the level of non-local operators. The part of the string operator corresponding to totally 
symmetric local tensor operators is projected out by 

[i>{X - z/2)-i,i,{X + z/2)Y''" = ^ / dti!{X-tz/2)z^{X + tz/2) (7) 







(we use the notation z = z'^^a). The subtraction of traces in the local operators implies that the twist-2 string 
operator contracted with z^ should satisfy the d'Alembert equation with respect to z: 



U,[i^{X-tz/2)z^{X + tz/2)\'"'"' ' = 0. (8) 
Transversality and twist-3 operators. In the coordinates X and z, the transversality conditions (^) are 

^n{,.}(z|x) = i^n[^.](^W ' ai;n[^H(^W = ^a|;;n{^^}(^m • (9) 

Consider the part of the current product given by Eq. (^) with the string operators replaced by their twist-2 parts. 
From Eq. (|) and {3 / dzp)[zp/ {2n^ z^)] = ~i5^^\z) it follows that [d / dzp)n^;^^f^ = 0, (c>/9z^) n{™f ^ q 
Since forward matrix elements are zero for all total derivative operators, this guarantees the transversal- 
ity of the twist-2 contribution in the case of deep inelastic scattering. In the non-forward case, we have 
(5/9X^)n*^'j*~2 ^ 0, (9/aX^)n*^^^*-2 7^ O , and (|) is violated. The non-transverse terms in the twist-2 contri- 
bution can only be compensated by contributions from operators of higher twist. In fact, the necessary operators 
are contained in the part of the string operator which was dropped in taking the twist-2 part. Incorporating 
QCD equations of motion, it is possible to show ||l^ that the twist> 2 part involves the total derivatives of string 
operators 

o /• 1 

^(-z/2)7„V(2/2) - [Vi(-z/2)7c<^(z/2)]''"" = ^eaCpKZCTTT^ / dtt i^{-tz /2)-i^-i5^j{tz/2) + . . . . (10) 
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The ellipses stand for quark-gluon operators (we do not write them explicitly since they are not needed to restore 
transversality of the twist-2 contribution). The relation for the operator with Dirac matrix 7075 is obtained by 
changing 7^ — > 7a 75, 7^75 — > 7k- The operators appearing under the total derivative on the R.H.S. of Eq. ( [l0| ) 
and its 7q75 analog are still the full string operators with no definite twist. Hence, one can decompose them into 
a symmetric (i.e., twist-2) part and total derivatives, and so on; thus expressing the original string operator as 
the sum of its symmetric part and an infinite series of arbitrary order total derivatives of symmetric operators. 
This series can be summed up in a closed form (the details of our calculation are presented elsewhere [ p5| ; similar 
expressions were derived independently in [p|jic|]). Up to operators whose matrix elements give 0{t) contributions 



to the Compton amplitude, the result is 

iv 
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/o Jv 







dXa 

i;{^vz/2)zj5^j{vz/2) + ... . (11) 



^13 • 

An analogous formula applies to the operators with 70- — > 70-75; one should just replace z —>■ £75, £75 z. 

III. PARAMETRIZATION OF NONFORWARD MATRIX ELEMENTS 

Double distributions. To get the amplitude for deeply virtual Compton scattering off a hadronic target we need 
parametrizations of the hadronic matrix elements of the uncontracted twist-2 string operators O^, O^a appearing 
in Eq. (^. We will derive them from Eq. (pd|). For simplicity, we consider here one quark flavor and the pion 
target, which has zero spin and practically vanishing mass. In this case, the matrix element of the contracted 
axial operator z'^O^a{z \ 0) (parametrized in the forward limit by the polarized parton density) is identically zero. 
Thus we need only the parametrization for the matrix element ip — r /2 \ 0{z | 0) | p + r/2) of the contracted vector 
operator 0{z \ 0) = z'^Oa{z \ 0). With respect to z, it can be regarded as a function of three invariants {pz), (rz) 
and z'^. For dimensional reasons, the dependence on is through the combinations tz'^ and p'^z'^ only. Since we 
are going to drop 0{t) and 0(j>^) terms in the Compton amplitude, we ignore the dependence on z'^ and treat 
this matrix element as a function of just two variables (pz) and (rz). Incorporating the spectral properties of 
nonforward matrix elements |l^ , we write the plane wave expansion in the form 

(p~r/2\O{z\0)\p + r/2) = 2{pz) f dif ^ ^-'^'''''^ f{i, a) da + (rz) f e^'"'^''''^''^ D{a) da , (12) 

where k = ip + ar/2, f{x,a) is the double distribution (DD) and D{a) is the Polyakov- Weiss (PW) distribution 
amplitude absorbing the (pz)-independent terms. From this parametrization, we can obtain the matrix elements 
of original uncontracted string operators, (pi]), including the kinematical twist-3 contributions. We consider first 
the part coming from the double distribution term in Eq. (|l2|); the contributions from the PW-term will be included 
separately. In matrix elements, the total derivative turns into the momentum transfer, id/dX„ r„ ~ 2^p^ + A^-. 
Similarly, we write k = (x + £,a)p + aA/2. This gives 

^{p-r/2\O^{z\0)\p + r/2) = f dif ' ' da /(i, a) (]5<,e''(*+«")(^'^)"'"('^^)/2 

2 J-i J-l+\j:\ I 

+ i [A„{pz) - p^{Az)] d^,„e~'''(*+«")(''^)-™"(^^)/2 [sin(5(rz)/2) -mcos(w(rz)/2)]| . (13) 

Skewed distributions. Expanding exp[— iQ;(Az) /2] = 1 — ia{Az)/2-\- . . . and keeping only terms up to those linear 
in the transverse momentuml^] A we get an expression in which the spectral parameter x appears in the exponential 



*Because of this truncation, the A^A^ terms in the expression for the amplitude T^i^ will be lost. If needed, they can be 
kept; see the discussion after Eq. (b2). 
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factors only in the combination x = x + ^a. Thus, we can introduce two skewed parton distributions: 

Hix,0 1 /■! /-i-l^l f 1 

> = / dx da6{x - X - ^a) f{x,a) < (14) 

A{x,^) J J-i J-i+\i\ [ a 

Note that, in our case, the DD /(i, a) is even in a and odd in x. As a resuh, the functions H and A satisfy the 
symmetry relations 

H{x,0 = -H{~x,0 , H{x,0=H{x,-0 , A{x,0=A{-x,0 , A{x,0 = -A{x,-0. (15) 

Furthermore, because of the antisymmetry of the combination af{x, a) with respect both to x and a we have 

1 

dxA{x,0 = 0. (16) 







Hence, the distribution A{x, ^) cannot be a positive-definite function on < a: < 1. 

Uniting the cosine and sine functions with the overall exponential factor, e~™^(P^), one gets vx±v£^ combinations. 
Using (^5|), one can arrange that only vx + would appear: 



i (p - r/2 I a (z I 0) b + r/2) = / dx e'^^f 



H{x,i)^'-^A(x,0 



+ -[A„{pz)^p^{Az)] j dx[H{x,0 - A{x,^)] j dvv cos[{vx + vOipz)] . (17) 
In a similar fashion, we get parametrization for the matrix element of the axial string operator (^ij) : 

^ (p - r/2 I Or,„{z \0)\p + r/2) = ^ e^^p^ z^App^ J dx [H{x, " M^, 0] j^^ dvv sm[{vx + vS,){pz)] . (18) 



Note that it is expressed in terms of the same skewed distributions H(x^ ^) and A{x^ ^) which, in turn, are determined 
by the original double distribution f{x,a), see Eq. (Q). 



IV. DVCS AMPLITUDE FOR PION TARGET 



DD-generated contribution. Substituting the parametrizations ( |l7|) and (18) into Eq. 
Fourier integral over the separation z we obtain the Compton amplitude 



and performing the 



T — 



+ 



1 

1 



\xJ^ 

-1 x-^ 



iO 



, dxR^x^i) I d. (-^^+3ep,)A. 1 r ,,^(,^^) r,, A.(..+ep.) 



where i?(x,^) is a new SPD describing the kinematical twist-3 contributions: 



i?(x,0 



_dH{x,0 dA{x.,0 



dx 



dx 



(19) 



(20) 



AU 



This result for the Compton amplitude contains the same tensor structures as those obtained in Refs. 
three terms in Eq. ( |l9|) are individually transverse up to terms of order t,p^. 

Singularities. The first term is the twist~2 part with the tensor structure corrected exactly as suggested by 
Guichon and Vanderhaeghen |^]. The integral over x exists if H(x,£,) is continuous at x = ^, which is the case for 
SPD's derived from the DD's which are less singular than 1/i^ for i = and are continuous otherwise (see jl^). In 
particular, continuous SPD's were obtained in model calculations of SPD's at a low scale in the instanton vacuum 
[P_8| . The second term contributes only to the helicity amplitude for a longitudinally polarized initial photon. The 
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parameter integral over v gives the function [ln(a; + ^ — iO) — ln(2^ — iO)]/{x — ^) which is regular at a; = ^ and 
has a logarithmic singularity at a; = — ^. The integral over x exists if R{x,^) is bounded at a; = — ^, which again is 
the case in the DD-based models described in Ref. [0 . The third term of Eq. ( p^ corresponds to the transverse 
polarization of the initial photon. In this case, one faces the integrand l/[v{x — ^) + iO] which produces dv/v 
divergence for the v-integval at the lower limit. One may hope to get a finite result only if the integral 

m . /' (21) 

vanishes. From the definition of the skewed distributions H{x, ^) and A{x,^) it follows that 

dA{x,0 ^ dH{x,i) 
dx 9^ 

Hence, one can substitute R{x, ^) by the combination dH{x, ^)/dx + dH{x, ^)/d£, similar to that used in [|lO| within 
the context of WW approximation (in Refs. full SPD's are implied while our H(x,^) does not include the 

PW term). By analogy with Ref. [foi, we integrate the dH{x,S^)/dx term by parts. This gives 



dxJ^^ , (22) 



i.e., the ^ derivative of the twist-2 contribution. In general, the latter has a nontrivial ^-dependent form determined 
by the shape of SPDs (see, however, the discussion of the PW contribution below). We conclude that the twist-3 
part of the tensor amplitude T^i, diverges in case of the transverse polarization of the initial photon. A similar 
observation has been recently made in Ref. [|lO|. However, it is easy to see that the relevant tensor structure 
^fiiliy + S.Pi') is just a truncated version of the exactly gauge invariant combination A^(72,y which has zero projection 
onto the polarization vector €2 of the final real photon: (£292) = 0. The structure A^q2jj is obtained if one uses 
the original full form of the DD parametrization (^3|) . It appears from the term with the exponential factor of the 
argument —i[v{x + S.a){pz) + vS^{pz) + {va + v){Az)/2\ which is obtained by combining the sine/cosine functions and 
the exponential factor in the second term of Eq. (^). In the Compton amplitude, it gives rise to a contribution 
in which the argument of the quark propagator is (7 + v{x + ^a)p + v£,p + {va + v)A/2. Since (Ag), (Ap) and A^ 



are negligible, the denominator factors in Eq. (19) remain unchanged. In numerators, representing {va + v)A/2 



as [1 — (1 — a)v]A/2, we observe that A^((7^ + ^p^) converts into the A^{q^ + S^p^ + Ai,/2) = Af^q2^ term plus 



a vAf^A^ type contribution corresponding to a new SPD built from the (1 — a) f{x,a) DD (cf. (14)). Due to 



the extra v factor, the w-integral for the latter contribution is finite. Hence, for the physical DVCS amplitude, 
we find no evidence for factorization breaking in the kinematical twist-3 contributions, both in their 1/ \/ —q^ and 
1/q^ terms. It is quite possible that factorization breaks down at the 1/q^ level, but one needs to analyze 0(2"^) 
suppressed terms (i.e., twist-4 contributions) to see if it really happens. 

WW-type representation. We can express our results in another form, introducing new skewed distributions 
related to R{x, £,) via an integral transformation similar to that used by Wandzura and Wilczek |l9j [] . This 
allows us to establish the equivalence between our approach and the WW-type approximation proposed in Ref. |^] . 
Treating the combination xv + v£^ in (|l9|) as a new variable we define 

Rw{x, = R{y, dy C 5{yv + v^-x)dv ^ 0{x > dy - 9{x < dy . (23) 



y~^ J-i y-S. 

In terms of this transform, the matrix element of the vector operator ( p^ ) can be expressed as 

Rw{x,£,) ~ Rw{~x,£,) 



-{p-r/2\O,iz\0)\p + r/2)= / dxe-^-'^P'Up„H{x,0-^PAMMx,0 



(24) 



^WW-type integrals of parton distributions have originally appeared within the ^-scaling formalism 
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Note that only the odd part of Rw{x,S,) contributes here. In case of the axial operator (|l 



(25) 



only the even part of Rwix, ^) appears. The part of the Compton amplitude ( [l9| ) containing R{x, ^) can be written 
in terms of this new function as 



1 



2iqp) y_i 



A^'{q•' + ^p'') (q^ + 3^p^)A'' 
X - ^ + jO ^ x + ^-iO 



Rw{x, ^) dx 



(26) 



The integrals with l/(a;±^=piO) converge only if the function Rw{x, ^) is continuous for x According to Eq. 

(p3[), Rw{x,^) is given by the integral of R{y, 0/{y~0 from x to 1 if x > ^ and from x to — 1 if a; < ^. Evidently, 
X = — ^ is not a special point in the integral transformation (23), hence the function R\y(x,S,) is continuous at 
X = — ^. However, it is extremely unlikely that the limiting values approached by Rw{x,^) for x = ^ from below 
and from above do coincide. Indeed, the difference of the two limits can be written as the principal value integral 
(compare with ||lO|) 



-1 y-^ 



dy 



(27) 



which can be converted into the ^-derivative of the real part of the twist-2 contribution. This means that the 
singularity, which we observed as a straight divergence of the dv/v integral, in this approach appears due to an 
unavoidable discontinuity of the Rwix,S,) transform at x = ^. 

Contribution from the PW-term. The contribution of the PW term to the vector operator 



1 



(p-r/2|a(0|z)|p + r/2) 



PW-Torm 



(28) 



has a simple structure corresponding to a parton picture in which the partons carry the fractions (1 ± q;)/2 of 
the momentum transfer r. Since only one momentum r is involved, this term can contribute only to the totally 
symmetric part of the vector string operator: it "decouples" in the reduction relations (p^). In particular, the 
PW term does not contribute to the second contribution in Eq. (|l|) which is generated by decomposition of the 
axial string operator: both derivatives, with respect to X and z, give rise to the momentum transfer r, whence the 
contraction with the e-tensor in (^ij) gives zero. Thus, the PW-contribution should be transverse by itself. Indeed, 
a straightforward calculation gives 



T, 



1 

{rq) 



rtiqu + qtiry - gi^u{rq) + r^r^ 



da , 



(29) 



which evidently satisfies g^T^i/lpw = 0, r^T^i/lpw = 0. Hence, this term can be treated as a separate contribution. 

Alternatively, one may include it into the basic SPD H{x, £,) and all SPD's derived from H{x, £,). Specifically, for 
^ > 0, the PW contribution to ff(x,0 is D{x/^)0{\x\ < £,) (l6); it contributes - x)D' {x/^) 9{\x\ < ^)/^2 [^here 
D'{a) = {d/da)D{a)] to i?(x,^); furthermore, the PW contribution to Rw{x,^) is r'(x/^)6'(|x| < 0/^. Inserting 
these functions into Eqs. ( |ig|) a nd ( p6|) one rederives Eq. (p9|). One can also observe that the PW term gives zero 
contribution into /(^), Eq. (^1[). 



V. CONCLUSIONS 



In this paper, we have studied the DVCS amplitude T^j, making use of the light-cone expansion in terms of QCD 
string operators in the coordinate space. We have demonstrated that transversality of the light -cone expansion can 
be maintained by including a minimal set of "kinematical" twist-3 operators, which appear as total derivatives of 
twist-2 operators. Incorporating the formalism of double distributions, we established that the kinematical twist-3 
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contributions are described by a skewed parton distribution Rw(x,£,) which can be derived from the basic twist-2 
double distribution /(i,a). The new SPD R\y{x,£,) has the structure of a generahzed WW transform. Due to 
discontinuities of Rwix,S,) ai x — ^, the factorization for the tensor amphtude T^i, breaks down at the twist-3 
level for the part corresponding to the transverse virtual photon. However, this term gives zero projection onto 
the polarization vector of the final real photon. 
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